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Abstract Unlike the Boussinesq, KdV and BBM equations, the celebrated Casamma- 
Holm ( CH) equation can model both phenomena of soliton interaction and wave break- 
ing. Especially, it has peaked solitary waves in case of u = 0. Besides, in case of 
u ^ 0, its solitary wave "becomes C°° and there is no derivative discontinuity at its 
peak", as mentioned by Camassa and Holm However, it is found in this article 
that the CH equation has peaked solitary waves even in case of u ^ 0. Especially, 
all of these peaked solitary waves have an unusual property: their phase speeds have 
nothing to do with the height of peakons or anti-peakons. Therefore, in contrast to the 
traditional view-points, the peaked solitary waves are a common property of the CH 
equation: in fact, all mainstream models of shallow water waves admit such kind of 
peaked solitary waves. 
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1 Introduction 

Since the solitary surface wave was discovered by John Scott Russell [16] in 1834, many 
models of solitary waves in shallow water have been developed, such as the Boussinesq 
equation [3], the Korteweg & de Vries (KdV) equation [9] and the Benjamin-Bona- 
Mahony (BBM) equation p], and so on. The Boussinesq and KdV equations are 
integrable and can model the soliton interaction of solitary waves and propagating 
waves with permanent form. However, they can not model breaking waves. The 
BBM equation has better analytic properties than the KdV equation, but it is not 
integrable, its traveling waves are not solitons and it can not model breaking waves 
(see [5]). In contrast to the KdV, Boussinesq and BBM equations, the celebrated 
Camassa-Holm (CH) equation [4] 

Ut + 2ujUx — Uxxt + 'iuUx = 2UxUxx + UUxxx (1) 

can model both phenomena of soliton interaction and wave breaking (see [5J), where 
u{x,t) denotes the wave elevation, x,t are the temporal and spatial variables, a; is a 
constant related to the critical shallow water wave speed, the subscript denotes the 
partial differentiation, respectively. Mathematically, the CH equation is integrable and 
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bi-Hamiltonian, thus possesses an infinite number of conservation laws in involution 
[1]. In addition, it is associated with the geodesic flow on the infinite dimensional 
Hilbert manifold of diffeomorphisms of line (see [5]). Thus, the CH equation ([T]) has 
many intriguing physical and mathematical properties. As pointed out by Fushssteiner 
[7], the CH equation ([T]) "has the potential to become the new master equation for 
shallow water wave theory" . 

Especially, when a; = 0, the CH equation ([1]) has the peaked solitary wave jl] 



When w 7^ 0, "the soliton solution of ([T]) becomes and there is no derivative 
discontinuity at its peak", as mentioned by Camassa and Holm [1]. As a result of it, 
investigations of the peaked solitary waves of the CH equation ([1]) were restricted in 
the case of u = (see [8] for example). 

Currently, the closed-form solutions of the peaked solitary waves of the KdV 
equation, the modified KdV equation, the Boussinesq equation and BBM equation 
are found by Liao [H]. In this article, it is found that the peaked solitary waves of 
the CH equation ([1]) also exist even in case of w 7^ 0. 



2 Peaked solitary waves of the CH equation in case 



Let us consider the propagating solitary waves of the CH equation with permanent 
form. Writing ^ = x — ct and w{C,) = c u{x,t), the original CH equation ([T]) becomes 



u{x, t) = cexp(— |x — ct\). 



of 7^ 




) 



w' + 3ww' = 2w'w" + WW' 



(2) 



subject to the boundary conditions 



w 0, 0, w" 0, as |^| +00, 



(3) 



where the prime denotes the differentiation. 
The linearized CH equation 




(4) 



has the peaked solitary wave 




provided co < c/2. 
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The solution of the equation ([2]) can be expressed in the form 

+00 

'iO = ^an exp(-n/i|^|), 



n=l 



where consant to be determined. Write A = w{0). Using the homotopy 

analysis method (HAM) pHHTB] , a analytic technique for highly nonlinear differential 
equations, it is easy to gain the series solution 

+00 

m=0 

Here 

u;o(0 = ^exp(-/i|e|) 
is the initial guess, Wm for m > 1 is governed by 

subject to the boundary condition 

Wm{0) =0, 0, as 1^1 +00, (8) 

where cq 7^ is an auxiliary parameter, called the convergence-control parameter, 

cf = r - f^'f 

is an auxiliary linear operator, and 

n 
j=0 

r 1 when n>l, 

^'^ \ otherwise. ^ ' 

The HAM is independent of small physical parameters. Especially, unlike other ana- 
lytic techniques, the HAM provides us a convenient way to guarantee the convergence 
of approximation series. For details, please refer to Liao [101 [121 113] • In fact, directly 
using the HAM-based mathematica package BVPh 1.0 (see Part II of [13]) for non- 
linear boundary- value/eigenvalue problems, it is straightforward to gain high-order 
analytic approximations of ([2|) and Q. For details, please refer to the Appendix. 

The accuracy of the mth-order approximation is defined by the averaged residual 
square of ([2|) in the domain ^ G [0, a]: 



a 



j=0 



d^, (11) 
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Table 1: m'(0+) of the mth-order analytic approximations and the corresponding 
residual squares of Camassa-Holm equation ([1]) in case of c = l,co = 1/4 given by the 
HAM-Based package BVPh 1.0 with the convergence-control parameter cg = —1. 



m 


A = 1/10 


A = 1/5 


«'(o+) 




n'(0+) 


'-'rn 


2 


-0.066733 


1.1 E-9 


-0.123744 


2.6E-7 


4 


-0.066668 


1.6 E-13 


-0.122575 


7.0E-10 


6 


-0.066667 


2.2 E-16 


-0.122484 


1.8E-11 


8 


-0.066667 


2.9 E-19 


-0.122475 


3.7E-13 


10 


-0.066667 


2.8 E-22 


-0.122475 


5.7E-15 



Table 2: u'{0+) of the mth-order analytic approximations and the corresponding 
residual squares of Camassa-Holm equation ([T]) in case ofc=l,a; = l/4 given by the 
HAM-Based package BVPh 1.0 with the convergence-control parameter cq = —1. 



m 


A = -1/10 


A = -1/5 


«'(o+) 


£ 


«'(o+) 


£ 

'-'rn 


2 


0.073804 


1.3 E-9 


0.152028 


3.4E-7 


4 


0.073854 


1.7 E-13 


0.152698 


7.4E-10 


6 


0.073855 


1.6 E-16 


0.152748 


9.3E-12 


8 


0.073855 


2.2 E-19 


0.152752 


2.1E-13 


10 


0.073855 


2.2 E-22 


0.152752 


3.4E-15 



where 

J\fw = w'" - fi^w' + 3ww' - 2w'w" - ww'". 

Since the wave elevation decays exponentially, we use a = 10 in this article. 

For simplicity, we study the case c = 1 henceforth. First, let us consider the case 
of a; = 1/4. As shown in Table [Tj the averaged residual squares of the lOth-order 
analytic approximations decreases to 2.8x10"^^ in case of A = 1/10 and to 5.7x10"^^ 
in case of A = 1/5, respectively. Besides, the corresponding values of m'(0+) (the 
limit is taken as — )■ from the right) quickly converge to -0.066667 and -0.122475, 
respectively. Similarly, as shown in Table |2l the averaged residual squares of the 
lOth-order analytic approximations decreases to 2.2x10^^^ in case of A = —1/10 and 
to 3.4x10"^^ in case of A = —1/5, respectively. These averaged residual squares are 
much smaller than those of numerical ones (see [8J as an example). Thus, without 
doubt, they are the solutions of the CH equation ([1]) in case of a; = 1/4. It is very 
interesting that all of the corresponding wave elevations are peaked solitary waves, 
as shown in Fig. [1] with peakons and Fig. [2] with anti-peakons for different values 
of A. Thus, the peaked solitary waves of the CH equation ([T]) indeed exists even 
in case of u ^ 0. Here, as pointed out by Constantin and Molinet [6], all of these 
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peaked solitary waves should be understood mathematically as weak solutions of the 
CH equation ([1]). However, physically, this kind of discontinuity of wave elevation 
widely appears in fluid mechanics, such as dam break [19] in hydrodynamics and 
shock waves in aerodynamics, which have clearing physical meanings. In fact, such 
kind of discontinuous problems belong to the so-called Riemann problem [2l [T5lfT8l fT9] , 
a classical field of fluid mechanics. 




X- c t 

Figure 1: The peaked solitary waves u{x,t) of the Camassa-Holm equation ([T]) in 
case of a; = 1/2 with the same phase speed c = 1. Black line: A = 1/10; Red line: 
A = 1/5; Blue line: A = 3/10; Green line: A = 2/5. 




Figure 2: The peaked solitary waves u{x, t) of the Camassa-Holm equation ([T]) in 
case of t<j = 1/2 with the same phase speed c = 1. Black line: A = —1/10; Red line: 
A = -1/5; Blue line: A = -3/10; Green line: A = -2/5. 

Secondly, let us consider the case of c = 1 and A = ±1/5 but different values 
of u. Similarly, using the HAM-based mathematica package BVPh 1.0 [13], it is 
straightforward to gain the corresponding convergent analytic approximations with 
high accuracy. It is found that all of them have peakons (when A = 1/5) or anti- 
peakons (when A = —1/5), as show in Fig. [3] and Fig. HI respectively. Note that, 
when a; = 0, the BVPh 1.0 gives the closed- form solution w = Aexp{—\x — ct\) of 
the peaked solitary wave, reported by Casamma and Holm [1]. It is found that, as u 
becomes larger, the wave elevation decays more slowly. This confirms once again that 
the CH equation ([T]) possesses the peaked solitary waves even when w 7^ 0. 



Figure 3: The peaked solitary waves u{x, t) of the Camassa-Holm equation ([T]) in 
case of A = 1/5 with the same phase speed c = 1. Red hne: u = 1/3; Green hne: 
Lo = 1/5; Blue line: u = 0. 




Figure 4: The peaked solitary waves u{x, t) of the Camassa-Holm equation ([T]) in 
case of y4 = — 1/5 with the same phase speed c = 1. Red line: u = 1/3; Green line: 
a; = 1/5; Blue line: u = 0. 
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3 Concluding remarks 

As mentioned by many authors, the CH equation ([1]) has many intriguing properties 
in physics and mathematics. Physically, it can model both phenomena of soliton 
interaction and wave breaking. Mathematically, it is integrable and bi-Hamiltonian, 
thus possesses an infinite number of conservation laws in involution and besides 
it is associated with the geodesic fiow on the infinite dimensional Hilbert manifold of 
diffeomorphisms of line. In addition, like Euler equation whose limiting wave also has 
a crest with a corner [T7], the CH equation admits peaked traveling waves not only for 
a; = but also for a; 7^ 0, as found in this letter. Especially, the phase speed of these 
peaked solitary waves has nothing to do with the height of the peakon or anti-peakon: 
this is an unusual property of the peaked solitary waves. 

Note that the closed-form solutions of peaked solitary waves of the KdV equa- 
tion, Boussinesq equation and BBM equation have been found by Liao [H]. Therefore, 
peaked waves seem to be common for the mainstream models of shallow water waves. 
Indeed, the elevation of these peaked solitary waves have discontinuity at crest. How- 
ever, this kind of discontinuity widely appears in fiuid mechanics such as dam break 
and shock waves which have clear physical meanings. In fact, such kind of discontin- 
uous problems belong to the classic Riemann problems. Therefore, the discontinuity 
of these peaked solitary waves are acceptable not only in mathematics but also in 
physics. It is true that such kind of peaked solitary waves have never been observed. 
So, it is an interesting and challenging work to observe them in laboratory. Certainly, 
further theoretical, numerical and experimental investigations on these peaked waves 
are necessary. 

Finally, it should be emphasized that, the discontinuity and/or singularity exist 
widely in natural phenomena, such as dam break in hydrodynamics, shock waves 
in aerodynamics, black holes in general relativity equation and so on. Indeed, the 
discontinuity and/or singularity are difficult to handle by traditional methods. But, 
the discontinuity and/or singularity can greatly enrich our understandings about the 
real world, and therefore should not be evaded. 
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Appendix 

The use of HAM-based Mathematica package BVPh 1.0 



Based on the HAM [9-11], the mathematica package BVPh 1.0 for nonlinear 
boundary- value/eigenvalue problems is developed and issued by Liao (Part II) [11], 
which is free available online. Using the BVPh 1.0, it is straightforward to gain the 
analytic approximations of the peaked solitary waves of the equations (2) and (3) for 
given w{0) = A. Here, we briefly describe how to do it in case of c = 1,0; = 1/4 and 
A = 1/10. 



1. First, download the BVPh 1.0 (the code flle is named by BVPh_1.0.txt) online 



( http://numericaltank.sjtu.edu.cn/BVPh.htm ) and save it in a directory such 



as C:/math/CH as an example. 

2. Then, run the computer algebra system Mathematica, and type the following 
command one by one: 



SetDirectory ["C:\math\CH"] ; 



«InputCH.txt 



The flle named InputCH.txt contains the following Mathematica commands and 
necessary deflnitions for BVPh 1.0: 



(* Install the BVPh 1.0 *) 
«BVPhl_0.txt; 



(* Define the physical and control parameters *) 
TypeEQ = 1 ; 

ApproxQ = ; 

ErrReq = 10" (-30); 

zRintegral = 10; 

(* Define the governing equation *) 
mu2 = l-2*omega/c ; 

f [z_,u_,lambda_] := D [u, {z ,3}] -niu2*D [u,z] \ 

+ 3*u*D[u,z] - 2*D[u,z]*D[u,{z,2}] - u*D[u,{z,3}] ; 



(* Define Boundary conditions *) 
zR = Infinity; 
OrderEQ = 3; 

BC[l,z_,u_,lambda_] := Limit[u-A, z -> ] ; 
BC[2,z_,u_,lambda_] := Limit [u, z -> zR ] ; 
BC[3,z_,u_,lambda_] := Limit [D [u,z] , z -> zR ] ; 



(* Define initial guess *) 

mu = Sqrt [mu2] ; 

u [0] = A*Exp [-inu*z] ; 

(* Define output term *) 

output [z_,u_,k_] := Print ["output = ",D[u[k],z] /. z->0//N] ; 

(* Defines the auxiliary linear operator *) 
L[u_] := D[u,{z,3}] - mu2 * D[u,z]; 

(* Print input and control parameters *) 
Print Input [u [z] ] ; 

(* Set convergence-control parameter cO and physical parameters 
cO = -1 ; 
A = 1/10; 
omega = 1/4; 
c =1; 

Print [" cO = ",cO, " omega = omega, " c = ", c, " A 

(* Gain up to lOth-order HAM approximation *) 
BVPh[l,10] ; 

(* Get results in the whole domain *) 
For[k=0,k<=10,k++,W[k] = U[k] /. z-> Abs[x]]; 

(* Show the 5th and lOth-order approximation *) 

Plot[{W[5] ,W[10]},{x,-10,10>,PlotRange->{Min[A,0] ,Max[A,0]}] 



